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Abstract — Several recent studies have found that clouds are sur- 
rounded by a transition zone of rapidly changing aerosol optical 
properties and particle size. Characterizing this transition zone 
is important for better understanding aerosol-cloud interactions 
and aerosol radiative effects, and also for improving satellite 
retrievals of aerosol properties. This letter presents a statistical 
analysis of a monthlong global data set of Cloud-Aerosol Lidar 
and Infrared Pathfinder Satellite Observation (CALIPSO) lidar 
observations over oceans. The results show that the transition zone 
is ubiquitous over all oceans and extends up to 15 km away from 
clouds. They also show that near-cloud enhancements in backscat- 
ter and particle size are strongest at low altitudes, slightly below 
the top of the nearest clouds. Also, the enhancements are similar 
near illuminated and shadowy cloud sides, which confirms that the 
asymmetry of Moderate Resolution Imaging Spectroradiometer 
reflectances found in an earlier study comes from 3-D radiative 
processes and not from differences in aerosol properties. Finally, 
the effects of CALIPSO aerosol detection and cloud identification 
uncertainties are discussed. The findings underline the importance 
of accounting for the transition zone to avoid potential biases in 
studies of satellite aerosol products, aerosol-cloud interactions, 
and aerosol direct radiative effects. 

Index Terms — Aerosol, cloud, Cloud-Aerosol Lidar and 
Infrared Pathfinder Satellite Observation (CALIPSO), lidar, 
satellite. 

I. Introduction 

S EVERAL recent studies have indicated that clouds are 
surrounded by a transition zone of rapidly changing aerosol 
optical properties and particle size (e.g., [l]-[5]). Characteriz- 
ing this transition zone is important for better understanding 
two critical yet poorly understood aspects of anthropogenic cli- 
mate change — aerosol-cloud interactions and aerosol radiative 
effects — [6], and also for devising effective sampling strategies 
for satellite measurements of aerosol properties. 

Satellites offer excellent opportunities to study the transition 
zone near clouds, although the separation between clouds and 
cloud-free areas is often uncertain, and ambiguities can occur 
at all scales [7]. As Charlson et al [8] summed up, “Multiple 
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lines of evidence exist that call into question the degree to 
which clear and cloudy sky can be separated. . Consequently, 
improving our knowledge of the transition zone requires im- 
provements in understanding and separating remote sensing 
uncertainties and physical phenomena such as aerosol swelling 
in humid air. 

One reason why such improvements are important is that 
aerosol changes and remote sensing uncertainties near clouds 
create a dilemma: Excluding the transition zone in order to 
avoid its remote sensing uncertainties can bias a study toward 
low aerosol optical depths, thus weaker radiative effects cal- 
culated from these optical depths, but including the transition 
zone despite the remote sensing uncertainties can bias the study 
toward too high aerosol optical depths and stronger radiative 
effects. 

Lidar instruments offer excellent opportunities for studying 
the transition zone, because their measurements are not affected 
by 3-D cloud radiative effects or instrument blurring (e.g., 
[9]— [12]). Recently, Tackett and Di Girolamo [13] pioneered 
the use of the spacebome Cloud-Aerosol Lidar and Infrared 
Pathfinder Satellite Observation (CALIPSO) lidar for studying 
the transition zone. Their study of aerosols over the Caribbean 
Sea showed that aerosol properties change rapidly near clouds 
both because aerosols swell in the humid environment and 
because the collision and coalescence of aerosol-containing 
cloud droplets and their subsequent evaporation produce fewer 
but larger aerosol particles. Additionally, aerosol nucleation by 
aqueous-phase reactions is also a factor. 

While the pioneering study of Tackett and Di Girolamo [13] 
focused on areas up to 3-4 km away from clouds in the Carib- 
bean, this letter analyzes a large global data set of CALIPSO 
measurements, which can reveal even the subtle changes that 
occur at the more distant portions of the transition zone. 

n. Data and Methodology 

This letter analyzes a monthlong (from September 15 to 
October 14, 2008) global data set of CALIPSO lidar data over 
all oceans free of sea ice. Our statistical analysis uses the 
operational 1-km-resolution CALIPSO cloud mask and 333-m- 
resolution attenuated backscatter profiles at both 532- and 
1064-nm wavelengths. Most results shown are for nighttime 
data simply because the weaker background illumination cre- 
ates smaller noise and allows more reliable cloud detection, but, 
as it will be illustrated in the last figure, daytime results showed 
very similar behaviors. 

Following Tackett and Di Girolamo [13], we characterize the 
transition zone by examining how median attenuated backscat- 
ter values vary with distance from the nearest cloud. Median 
values are well suited for this purpose because they are less 
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Fig. 1. (a) Global median vertically integrated attenuated backscatter for all nighttime CA1JPSO observations over ice-free oceans. (Inset) Cumulative histogram 

of cloud-free profiles, (b) Median color ratio. These and subsequent plots are based on a monthlong data set in September-October 2008. The distances in all 
figures are to the nearest cloud below 3-km altitude. The error bars in this and subsequent figures indicate uncertainties estimated using the bootstrap method [19]. 
Because this method assumes that all profiles are independent samples, it represents a lower limit to actual uncertainties. 


affected than mean values by uncertainties in cloud detection 
and by unobserved clouds that lie right near the satellite track. 
We calculate median values by identifying the 50th-percentile 
points in high-resolution backscatter and color ratio histograms. 
Let us note that, although this letter shows only median values, 
mean values show generally similar but noisier behaviors. 

This letter focuses on clear areas closest to low-level clouds 
(below 3 km), but we note that results for clear areas near any 
clouds were very similar to those presented here. One reason 
for the similarity appears in the inset of Fig. 1(a), which shows 
the cumulative histogram of clear-sky profiles as a function of 
distance to the nearest cloud that has a cloud top below 3 km. 
The figure shows that 50% of all clear profiles are closer than 
4.4 km to a cloud below 3 km, and 75% of them are closer than 
20 km. The remaining 25% of clear profiles include both those 
clear profiles that are closest to a cloud higher than 3 km and 
those clear profiles that are more than 20 km away from any 
cloud. We note that Twohy et al [5] found 66% of all clear 
areas they observed over the Indian Ocean to be closer than 
4 km to clouds. 

Our results indicate that much more clear profiles lie close 
to clouds below 3 km than to higher clouds, even though 
the highest cloud top is more often above than below 3 km 
(40% versus 38% of all CALIPSO profiles in our data set). 
This occurs because high clouds often form large continuous 
layers, and so, only a few clear profiles can be close to high 
clouds (at the two ends of the large cloud decks) — whereas low 
cloud layers often consist of numerous small clouds separated 
by clear gaps, which allows many clear profiles to be close 
to them. The CALIPSO cloud detection algorithm [14] may 
further emphasize this tendency when it reduces the effects of 
observational noise by detecting faint (and often high) clouds 
after averaging observed backscatter values over larger scales, 
which fills in any small gaps in thin high clouds. 

Finally, we note that we show results of attenuated backscat- 
ter values integrated up to 3-km altitude and that the results 
behaved very similarly when integrating to 15-km altitude. We 
chose to present the results for low-level backscatter because 
this is where the bulk of near-cloud aerosol changes occur. 
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Fig. 2. Map of differences between median vertically integrated 532-nm 
backscatter values ( p ) at different distances from the nearest cloud. 
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Fig. 3. Altitude dependence of median enhancements near clouds, (a) 532-nm attenuated backscatter versus altitude above sea level, (b) Color ratio versus 
altitude above sea level, (c) 532-nm attenuated backscatter versus altitude above the top of the nearest cloud, (d) Color ratio versus altitude above the top of the 
nearest cloud. 


m. Results on Particles Near Clouds 

Fig. 1(a) confirms the overall tendency found in earlier pa- 
pers that attenuated backscatter increases systematically toward 
clouds. The figure indicates that this tendency is sufficiently 
ubiquitous to appear in global statistics and that the transition 
zone typically extends to about 15 km. 

Fig. 1(b) shows that the spectral behavior of attenuated 
backscatter (/?) also changes near clouds: The median “color 
ratio” (which is slightly smaller than the ratio of medians at 
the two wavelengths but follows the same behavior) increases 
systematically near clouds. Because backscatter decreases with 
wavelength faster for smaller particles than for larger particles 
(i.e., smaller particles have larger Angstrom exponents), the 
color ratio increase in Fig. 1(b) indicates that particle size 
increases near clouds. This increase is consistent with the 
effects of both aerosol humidification and cloud detrainment, 
but establishing the relative importance of these two processes 
as in [13] will require a more detailed study. 

Fig. 2 shows that the increase near clouds in Fig. 1 is a 
global phenomenon and that the overall trend is neither caused 
by a single dominant region nor is it an artifact of combining 
observations from vastly different regions of the Earth. While 
undoubtedly there are regional variations, much of the varia- 
tions in Fig. 2 likely come from sampling issues. 


Fig. 3 shows the vertical structure of near-cloud backscat- 
ter enhancements by comparing the altitude dependence of 
backscatter at four intervals of distance to clouds (0-5, 5-10, 
10-15, and 15-20 km). Panels (a) and (b) of the figure show 
that aerosol content and size (i.e., 532-nm backscatter and 
color ratio) both decrease with altitude and that the near-cloud 
enhancements also decrease with altitude. The exception is a 
small increase at very low altitudes, a trend most pronounced at 
tropical areas. 

Panels (c) and (d) of the figure show the way backscatter 
and color ratio depend on altitude relative to the top of the 
nearest cloud. They reveal that aerosol content and size drop 
sharply above cloud altitude, as convection is not as effective 
in moving aerosols from low levels to above cloud tops. Near- 
cloud enhancements are also much larger below than above 
the cloud top. The maximum enhancement at cloud altitude is 
consistent with cloud detrainment occurring here by default, 
with aerosol humidification being strongest when humidity 
approaches saturation (e.g., [5]). 

IV. Effects of Remote Sensing 
Limitations Near Clouds 

While aerosol particles are responsible for a large part of 
the near-cloud enhancements discussed earlier, misidentified or 
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Fig. 4. Results on remote sensing uncertainties near clouds, (a) and (b) Median vertically integrated 532-nm attenuated backscatter and color ratio for nighttime 
CALIPSO observations over ice-free oceans. These panels consider separately the atmospheric columns that are free of both cloud and aerosol layers and the 
columns that contain an aerosol layer, as specified by the operational CALIPSO product. The inset in (b) shows the vertical profile of median color ratio in 
atmospheric columns that are 15-20 km away from the nearest detected clouds and where the operational CALIPSO product did not detect any aerosol layers, 
(c) and (d) Median near-cloud vertically integrated 532-nm attenuated backscatter and color ratio in daytime observations over all oceans between 30° and 60° 
latitudes in both hemispheres. The solid and dashed lines are for profiles that are on the sunlit and shadowy sides of the closest cloud, respectively. 


undetected cloud particles are also likely to contribute. Unde- 
tected cloud particles contribute to clear-sky statistics because 
the separation of cloudy and cloud-free areas is sometimes 
ambiguous due to uncertainties both in detecting layers of 
scattering particles and in distinguishing clouds from aerosols 
[14]. Identifying cloud particles is particularly difficult when 
small clouds cover only a portion of the lidar field of view — a 
situation that often occurs in sparse cumulus fields, where 
Koren et ah [7] found very small clouds to be fairly abun- 
dant. Although using medians instead of mean values greatly 
reduces the influence of cloud detection errors, we note that 
an ideal cloud mask would slightly weaken the near-cloud 
enhancements. Another aspect of cloud detection uncertainties 
is that clouds that lie next to the lidar path do not get detected, 
which causes an overestimation of the distance to the nearest 
cloud. This results in somewhat overestimating the width of the 
transition zone. 

As an example of remote sensing issues. Fig. 4(a) and (b) 
shows that lidar backscatter and color ratio increase near clouds 
even in pristine atmospheric columns where the operational 
CALIPSO product did not detect any aerosol layers. In com- 
parison to all cloud-free profiles, the increase in pristine areas 


is particularly well pronounced within 4-5 km from the clouds 
identified by CALIPSO. An unexpected feature of Fig. 4(a) 
is that backscatter right next to clouds is higher for pristine 
areas than for areas with detected aerosol layers. A possible 
reason why aerosol is not detected despite the higher backscat- 
ter values is that, in the coarse-resolution (>5-km) CALIPSO 
aerosol detections, lower “pristine area” backscatter farther 
away from clouds pushes area-averaged backscatter values be- 
low the aerosol detection threshold. The inset in Fig. 4(b) shows 
that color ratios below the cloud top greatly exceed the value of 
~0.07 that is expected from clear air molecular scattering (e.g., 
[14]). These results indicate that aerosol or cloud particles are 
often present even in these pristine profiles. 

Therefore, while particle backscatter and size certainly in- 
crease near the clouds detected by CALIPSO, both the nature of 
these particles (cloud versus aerosol) and the processes creating 
them (e.g., aerosol swelling in humid air) need to be clarified 
using a quantitative assessment of remote sensing limitations in 
particle detection and classification. 

Finally, Fig. 4(c) and (d) confirms the importance of remote 
sensing effects that influence satellite measurements of solar 
reflectance near clouds. In particular, the figures show daytime 
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near-cloud enhancements of lidar backscatter between 30° and 
60° latitudes. Each panel contains two curves: one for clear- 
sky profiles that are on the sunlit side of the nearest cloud 
(toward the Equator from the cloud) and the other that are on 
the shadowy side of the nearest cloud (toward the poles from 
the cloud). The figure indicates that the increase is very similar 
on both sides of clouds. This confirms that the increase being 
larger near sunlit cloud sides in Moderate Resolution Imaging 
Spectroradiometer (MODIS) solar reflectance observations is 
most likely a result of 3-D radiative processes [12], as opposed 
to asymmetric changes in aerosol properties. This implies that 
3-D radiative effects are important in shaping clear- sky radia- 
tion fields near clouds and should be accounted for in interpret- 
ing passive space-based aerosol measurements. 

V. Conclusive Thoughts 

This letter presents a statistical analysis of a monthlong 
global data set of CALIPSO lidar observations over oceans. The 
results show the following. 

• The transition zone of enhanced backscatter and particle 
size around clouds is ubiquitous over all oceans. 

• The transition zone is much wider (~15 km) than the 
ranges examined in most earlier studies. 

• Near-cloud enhancement is strongest at low altitudes, 
slightly below the top of the nearest clouds. 

In addition, the results also provide information on uncertain- 
ties in aerosol remote sensing near clouds. 

• CALIPSO enhancements are similar near illuminated and 
shadowy cloud sides, confirming that the differences found 
in an earlier study of MODIS solar reflectances [12] 
indeed come from 3-D radiative effects and not from 
differences in aerosol properties. 

• Aerosol or cloud particles cause backscatter and color ratio 
enhancements near clouds even in pristine profiles where 
the CALIPSO products did not detect any aerosol or cloud. 

The width and ubiquity of the transition zone near clouds 
imply that studies of aerosol-cloud interactions and aerosol 
direct radiative effects need to quantitatively account for aerosol 
changes near clouds. These changes can cause systematic bi- 
ases toward smaller aerosol optical depths, thus weaker ra- 
diative forcings calculated from these optical depths, if areas 
near clouds are not included in order to avoid aerosol retrieval 
ambiguities that arise near clouds (e.g., 3-D radiative effects, in- 
strument blurring, and cloud contamination). On the other hand, 
including aerosol products near clouds despite these ambigui- 
ties overestimates aerosol radiative forcing. As a result, there 
is an urgent need for developing methods that can assess and 
account for remote sensing challenges (e.g., for passive remote 
sensing, [15], [16]) and thus allow for including the transition 
zone into the study. Alternatively, one may remove biases by 
seeking a balanced sampling strategy that may compensate for 
excluding both near-cloud and far-from-cloud areas to a similar 
degree (e.g., [17]). Finally, statistical adjustments (e.g., [18]) of 
results limited to areas far from clouds may also be possible 
if the statistical behavior of the transition zones becomes well 
understood. More detailed future analysis of transition zones 
is expected to provide further insights into their origins and 
influence on the atmosphere, radiative processes, and satellite 
data products. 
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